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High energy emission from AGN cocoons in clusters of galaxies 
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Gamma-ray emission from cocoons of young radio galaxies is predicted. Considering the process of adiabatic injection of 
the shock dissipation energy and mass of the relativistic jet into the cocoon, we find that the thermal electron temperature 
of the cocoon is typically predicted to be of the order of ~ MeV, and is determined only by the bulk Lorentz factor of the 
jet. Together with the time-dependent dynamics of the cocoon expansion, we find that young cocoons can yield thermal 
Bremsstrahlung emissions at energies ~MeV. Hotter cocoons (i.e., GeV) for younger sources are also discussed. 
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1 Introduction 

Relativistic jets in active galactic nuclei (AGNs) are widely 
believed to be the dissipation of kinetic energy of relativistic 
motion with a Lorentz factor of order ~ 10 produced at the 
vicinity of a super-massive black hole at the galactic cen- 
ter (Begelman, Blandford and Rees 1984 for reviews). The 
jet in powerful radio loud AGNs (i.e., FR II radio sources) 
is slowed down via strong terminal shocks which are iden- 
tified as hot spots. The shocked plasma then expand side- 
ways and envelope the whole jet system and this is so called 
a cocoon or a bubble (Fig. 1). The cocoon is a by-product 
of the interaction between AGN jets and surrounding intra- 
cluster medium (ICM). The internal energy of the shocked 
plasma continuously inflates this cocoon. So far little atten- 
tion has been paid to observational feature of the cocoon, 
since they are usually invisible in GHz bands because of the 
synchrotron cooling for older electrons. As a result, we 
just see a part of the cocoon. The visible part is so-called ra- 
dio lobes in which relatively fresh electrons are filled in. In 
Fig. 2, we show one good sample of cocoon emission from 
the powerful radio galaxy Cygnus A for the evidence of its 
existence. 



* Corresponding author: e-mail: motoki.kino@nao.ac.jp 

1 On the contrary, the emission from the shell made of the shocked ICM 
(see Fig. 1) has been explored by many authors (e.g., Heinz, Reynolds and 
Begelman 1998; Sutherland and Bicknell 2007). Since the shells have non- 
relativistic velocities, the emission from the shells is predicted in the X-ray 
band. 

2 The data obtained by Carilli et al. (1991) was re-analyzed to obtain the 
map shown in Fig. 2. The observation was carried out with VLA A config- 
uration at 330 MHz on 1987 August 18. The data analyzes is performed by 
standard manner with Astronomical Image Processing System (AIPS). The 
flux-scale is determined by comparison with 3C 286 and 3C 48 using the 
AIPS task SETJY. The image is obtained with the DIFMAP after a number 
of self-calibration iterations. 



Among a variety of AGN bubbles, a population so called 
compact symmetric objects (CSOs) has been widely inves- 
tigated in various ways (e.g., Fanti et al. 1995; Readhead 
et al. 1996; O'Dea & Baum 1997; de Vries et al. 1997; 
O'Dea and Baum 1998; Stanghellini et al. 1998; Snellen et 
al. 2000; Dallacasa et al. 2000; Giroletti et al. 2003; Nagai 
et al. 2006; Orienti et al. 2007; Kawakatu et al. 2008). CSOs 
are smaller than 1 kpc and the previous studies support the 
youth scenario in which CSOs propagate from pc scales 
thrusting away an ambient medium and growing up to FR 
II radio galaxies. CSOs are thus recognized as newly born 
AGN jets, and they are crucial sources to explore physics of 
AGN bubbles in their early days. 

In this study, we propose that "young AGN bubbles" 
are a new population of 7-ray emitters in the Universe. The 
layout of the paper is as follows. We review the expand- 
ing cocoon model in §2 following our previous works (Kino 
and Kawakatu 2005; Kino, Kawakatu and Ito 2007, KKI07 
hereafter) We show the predicted MeV gamma emission 
from young cocoons in §3 . In §4, we further predict GeV 
gamma emission for smaller CSOs. Summary and discus- 
sion is given in §5. 

2 Cocoon inflation by exhausted jet 

Here we consider the time-evolution of an expanding co- 
coon inflated by the dissipation energy of the relativistic jet 
via terminal shocks. The adiabatic energy injection into the 
cocoon is assumed. Mass and energy conservation from the 
jet into the cocoon, which govern the cocoon pressure P c 
and mass density p c are written as 

% Pc(t)V c (t) _ 
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Fig. 1 A cartoon representation of interaction of the ICM 
with declining atmosphere and the relativistic jet in FR II 
radio galaxy. As a result, most of the kinetic energy of jet 
is deposited in the cocoon and it is inflated by its internal 
energy. 




Right Ascension (mas) 

Fig. 2 VLA image of Cygnus A at 330 MHz 
(A-configuration). The contour level starts with 
0.565 Jy beam -1 and increases from there by factors 
of 2. The convolving beam is 4.73 x 3.91 arcsec at the 
position angle of 10.1°. In this frequency, we can see the 
synchrotron emission from slightly colder electrons than 
those emitting GHz ranges. Hence we can identify the 
cocoon profile which are clearly different from familiar 
"classical double lobe". 



where %, V c , T® 1 , Jj and Aj, are the adiabatic index of the 
plasma in the cocoon, the volume of the cocoon, the kinetic 
energy and mass flux of the jet, and the cross-sectional area 
of the jet, respectively. The total kinetic energy and mass 
flux of the jet are Tj 01 = pjC 2 T 2 Vj, Jj = PjTjVj where p } , 
and Tj are mass density and bulk Lorentz factor of the jet 
(Blandford and Rees 1974). Hereafter we set Vj — c. The 
total kinetic power of the relativistic jet is defined as Lj = 
2Tj >1 (t)Aj (t) and it is assumed to be constant in time. 

As for the mass and kinetic energy flux of powerful rela- 
tivistic jets, numerical simulations tell us that no significant 



entrainment of the environmental matter takes place during 
the jet propagation (e.g., Mizuta et al. 2004). According to 
this, the mass and kinetic energy flux of the jet are regarded 
as constant in time. Then, the conditions of TJ 01 = const, 
and Jj = const leads to the relations of pj (t)Aj (t) = const 
and Tj (t) — const. In order to evaluate Lj, we use the shock 
jump condition of T?/9j = /3^ s picu (Kawakatu and Kino 
2006) where /3i ls (= «hs/ c ) an d Picm is the advance speed 
of the hot spot (3^ = 10~ 2 /3_ 2 an d the mass density of 
ICM, respectively. Using, the jump condition, Lj is given 
by 



Lj = 2 x 10 4 



Rt 



*-k pc 0-2 n -2 erg s (3) 
where we use Aj(i) = irR^ s (t), and the hot spot radius 
Rhs is given by i?k pc = i?hs(10 7 yr)/l kpc. As a fiducial 
case, we set the number density of the surrounding ICM as 
nicM(d) = p\cu{d)/m p = 10~ 2 cm _3 n_ 2 (<i/30 kpc) -2 
where d is the distance from the center of ICM and cocoon 
(see Fig. 1). Since the change of the index from —2 does 
not change the essential physics discussed in this work, we 
focus on this case for simplicity. Since Lj is the ultimate 
source of the phenomena associated with the cocoon, all of 
the emission powers which will appear in §3 should be less 
than Lj . 

The number density of total electrons in the cocoon is 
governed by the cocoon geometry and its plasma content. 
For convenience, we define the ratio of "the volume swept 
by the unshocked relativistic jet" to "the volume of the co- 
coon" as A(t). We denote V c (t) = 2(ir/3)Tl 2 Zl s (t), Z hs 
satisfies Zy^it") = /3hsc£, R c , and TZ = R c /Z^ s < 1 as the 
cocoon volume, the distance from the central engine to the 
hot spot, is the radius of the cocoon body, and the aspect- 
ratio of the cocoon, respectively (e.g., Kino and Kawakatu 
2005). Postulating that 1Z and Z\^j i?h s are constant in time, 



_ 2A i (t)v i t 



is evaluated as 



— Vc (t) 

A{t)^0AK- 2 Rl vc Z^pz\ (4) 

where Z30 = Zh s (10 7 yr)/30 kpc. Note that, in the case, 
the time dependence of A is deleted since V c oc t~ 3 and 
Aj oc i 2 . This case satisfies t>hs = const (e.g., Conway 
2002). The cocoon mass density p c (t) is controlled by the 
mass injection by the jet and it can be expressed as p c (t) w 
rjpj(£).4 = f^TT 1 p lGM (Z\y S {t))A where we use the shock| 
condition of T 2 pj = f3^ s picM- Adopting typical quanti- 
ties of FR II sources (e.g., Begelman, Blandford and Rees 
1984), the number density of the total electrons in the co- 
coon is given by 

-2 



5 (t) « 4 x 10- 5 An_ 2 r 10 /3 2 



t 



'(5) 



10 7 yr, 

where T = lOTio, and A = A/OA. Here we assume that 
the mass density of the pair plasma is heavier than that 
of electron-proton one, and then we adopt p c rs 2m e n e in 
the light of previous works (Reynolds et al. 1996; Wardle 
et al. 1998; Sikora and Madejski 2000; Kino and Takahara 
2004). The upper limit of thermal n e can be basically con- 
strained by the analysis of Faraday depolarization (Dreher 
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et al. 1987). However, the strong Faraday depolarization ob- ~ 10 3 times larger than that vF v ~ 10" 
served in CSOs (Cotton et al. 2003) are likely to be caused This is still less than the threshold of INTEGRAL. This 
by dense foreground matter such as narrow line region. There-|may be the reason for the lack of detection of MeV emis- 



fore n e in radio lobes of CSOs has not been clearly con- 
strained. 

Let us estimate the electron (and positron) temperature 
(T e ) and proton temperature (T p ). From Eqs. (HJ and d2j 
together with the equation of state P c ps 2n e kT e , we can 
directly derive the temperatures as 

kT e ps 1 Tio MeV, kT p ps 2 r 10 GeV (6) 

where we adopt the two temperatures condition of kT e ps 
(m e /m p )kT p . It should be stressed that the temperatures 
are governed only by Tj. It is also worth noting that the geo- 
metrical factors in Eqs. ((TJ and (O are completely cancelled 
out. One can naturally understand these properties by com- 
paring the well-established properties such as supernovae 
and GRBs. Constant temperature in AGN jet can be real- 
ized by the "continuous" energy injection into the expand- 
ing cocoon whilst temperatures of astrophysical explosive 
sources such as gamma-ray bursts and supernovae would be 
decreased because of "impulsive" injection of the energy. 
Thus the resultant temperatures are uniquely governed by 
Tj and they remain to be constant in time. 

3 MeV 7 emissions from a young cocoon 

3.1 Thermal MeV Bremsstrahlung emission 

The time-dependence of the thermal Bremsstrahlung lumi- 
nosity L Blcm is given by L B rom(^) °c n 2 e (t)Te /2 V c {t) oc 
tr 1 based on the cocoon expansion shown in the previous 
section. Hence it is clear that a younger cocoons are brighter 
Bremsstrahlung emitters than older cocoons. In a similar 
way, brighter synchrotron luminosity has been expected for 
younger radio galaxies (e.g., Readhead et al. 1996). With 
relativistic thermal Bremsstrahlung emissivity (Rybicki and 
Lightman 1979), the luminosity of the optically thin ther- 
mal Bremsstrahlung emission vL u at energies ~ 1 MeV is 
estimated as 

-l 



sion from young cocoons up to now. Fig. 3 shows that the 
XMM/Newton satellites can detect the low energy part of 
the thermal Bremsstrahlung from young cocoons. In MeV 
energy band, a proposed mission of detector SGD on board 
the NeXT satellite with the eye up to ~ 0.6 MeV (Taka- 
hashi et al. 2004) could detect the thermal MeV emission 
from those located slightly closer or younger with smaller 
Lorentz factor. 



3.2 Candidate sources 

In order to explore the extended cocoon emission in the 
X-ray band, one may think it is hard to distinguish over- 
lapping emission from the compact core of the AGN with 
limited spacial angular resolution of the current X-ray satel- 
lites. However, the averaged spectral index in X-ray band 
(T^) from the compact core of AGNs is softer than the 
Bremsstrahlung emission (Koratkar and Blaes 1999). Hence 
it is possible find candidate sources of the MeV cocoon by 
the value of Tx ■ As far as we know, there are two possi- 
ble candidates for the Bremsstrahlung emission. Those are 
B 1358+624 (Vinkat al. 2006) and PKS B 1345+ 125 (Siemigi-| 
nowska et al. 2008) actually shows Tx ~ L both of them 
were observed by XMM/Newton. Time variability of ob- 
served spectra is also the key to distinguish them. It is obvi- 
ous that the cocoon emission is constant in time whilst vari- 
ous emissions from the core of AGN should be highly vari- 
able. Hence steady emissions are convincingly originated in 
cocoons. 



iBrcm(t) « 2 X 10 4 ° n^Q 3 ^ 



10 



t 



10 7 yr 

Eq. ^ explains the reasons for the non-detection of the 
thermal emission from older cocoons. One is simply be- 
cause it is not very bright. The other is because the predicted 
energy range is ~ 1 MeV, the MeV-7 astronomy is still im- 
mature and it is sometimes called as "sensitivity gap" com- 
pared with the energy range below 10 keV and above GeV 
ranges (Takahashi et al. 2004). 

In Fig. 3, we show the predicted values of vF v for the 
cocoons with t = 10 7 yr and t = 10 4 yr located at the 
distance of D = 10 2 Mpc. The cocoon with t = 10 7 yr 
have vF v ~ 10 -14 ergcm _2 s _1 . The detection thresh- 
old of SPI instrument on board the INTEGRAL satellite 
is about vF v ~ 10 -9 erg cm _2 s _1 at ~ 1 MeV. For a 
young cocoon with t = 10 4 yr, the predicted luminosity is 



Intriguingly, "the diffuse X-ray emission" has been in- 
deed detected in PKS B 1345+ 125 with the size of the ex- 
tended emission is of order ~ 20 kpc by Siemiginowska et 
al. (2008). The diffuse emission might be associated with 
the radio lobes of this source. The X-ray emission is elon- 
gated towards the South- West similarly to the VLBI jet axis 
erg s _1 (7X e P ortec l by Stanghellini et al. (2001). If the emission is as- 
sociated with the radio emitting plasma, there are two pos- 
sibilities to explain the emission. Here we newly stress that 
the tail of the Bremsstrahlung emission could explain the 
emission. Based on Eq. ©, the X-ray luminosity of the 



emission shows Lx ps 1 X 10 erg s i . The observed Lx 



can be explained with the n e 



1 x 10 cm . Based 



on Eq. ©, the n e can be realized with [3h s ~ 10' 1 and 
niCM ~ 10 _1 cm' 3 for instance. Since the required n e is 
considerably large, the analysis of Faraday depolarization 
will be crucial for checking the upper limit of n e . The other 
possibility is non-thermal emissions from the lobes which 
is recently investigated by Stawarz et al. (2008). Observa- 
tions near the peak of thermal emission NeXT satellite will 
be crucial to distinguish whether the emission is thermal or 
non thermal one. 
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4 Younger radio sources as "hotter" bubbles 

So far, we discuss the cocoon property in the phase of no 
significant cooling. The phase roughly corresponds to medium| 
size symmetric objects (MSOs) and FR II galaxies. Since 
cooling timescales become shorter for smaller sources, cool-| 
ing effects for CSOs are more effective than the case for 
larger ones such as MSOs and FR lis. We consistently solve 
a set of equations which describes young bubble expansions 
including the effects of Bremsstrahlung emission and adia- 
batic loss together with the initial conditions indicated by 
CSO observations. Then we find that the bubbles have elec- 
tron temperature of ~ GeV at initial phases, the bubbles 
then cooled down to MeV by the adiabatic loss. We further 
estimate these 7-ray emissions and show that it could be 
detected with Fermi (GLAST) (Kino et al. 2009). 

5 Summary and discussion 

We have investigated the luminosity evolutions of AGN co- 
coons together with the dynamical evolution of expanding 
cocoon. Below we summarize the main results of the present 
work. 

1 . We newly predict the Bremsstrahlung emission peaked 
at MeV-7 band as a result of standard shock dissipation 
of relativistic jets in AGNs. The temperature of the co- 
coon is governed only by the bulk Lorentz factor of the 
jet Tj . The electron temperature T e relevant to observed 
emissions is typically predicted in the range of MeV for 
Tj ~ 10. Constant temperatures of plasma in the co- 
coon can be realized because of the continuous energy 
injection by the jet with constant Tj (KKI07). 

2. We further investigate younger bubble expansions in- 
cluding the effects of Bremsstrahlung emission and adi- 
abatic loss together with the initial conditions of CSOs. 
Then we find that the lobes initially have electron tem- 
perature of GeV and the lobes then cool down to MeV 
by the adiabatic loss. The 7-ray emissions could be de- 
tected with Fermi (GLAST) (Kino et al. 2009). 
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